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Warm dust in the terrestrial planet zone of a sun-like Pleiad: 
collisions between planetary embryos? 

Joseph H. Rhee^, Inseok Song^, B. Zuckerman^'^ 
ABSTRACT 

Only a few solar-type main sequence stars are known to be orbited by warm 
dust particles; the most extreme is the GO field star BD+20 307 that emits 
~4% of its energy at mid-infrared wavelengths. We report the identification of 
a similarly dusty star HD 23514, an F6-type member of the Pleiades cluster. A 
strong mid-IR silicate emission feature indicates the presence of small warm dust 
particles, but with the primary fiux density peak at the non-standard wavelength 
of ~9 nm. The existence of so much dust within an AU or so of these stars is not 
easily accounted for given the very brief lifetime in orbit of small particles. The 
apparent absence of very hot (>1000 K) dust at both stars suggests the possible 
presence of a planet closer to the stars than the dust. The observed frequency 
of the BD-l-20 307/HD 23514 phenomenon indicates that the mass equivalent of 
Earth's Moon must be converted, via collisions of massive bodies, to tiny dust 
particles that find their way to the terrestrial planet zone during the first few 
hundred million years of the life of many (most?) sun-like stars. Identification 
of these two dusty systems among youthful nearby solar-type stars suggests that 
terrestrial planet formation is common. 

Subject headings: circumstellar matter — infrared: stars — planetary systems: 
formation — stars: individual (HD 23514) — open clusters and associations: 
individual (Pleiades) 
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Introduction 



The Spitzer Space Telescope is now providing a wealth of new information about dusty 
stars in the Milky Way. However, because Spitzer is a pointed telescope, during its cryogenic 
lifetime it will examine only a modest portion of the sky, about 1%. For the very rare, very 
infrared-bright, nearby field star, the less sensitive Infrared Astronomical Satellite {IRAS) 
was actually the better search engine because it was an all-sky survey. Conse quently, we 



have an ongoing progra m to correlate the I RAS catalog with the Hipparcos (ISong et al. 



2OO2I : iRhee et al.ll2007bl ) and Tycho catalogs (IMelis et al.ll2008l ). To date we have identified 
at least three nearby stars of age ^100 Myr that emit at le ast a few percent of their energy 
at infrared wavelengths. These are: field star BD-l-20 307 (ISong et al.ll2005l: A. W einberger 
et al 2007, in preparation); a member of a field Tycho binary star (IMelis et al.ll2008l ): and the 
Pleiad HD 23514 (= HI 1132) that is the principal focus of the present paper. In addition 
to these three stars, BP Psc, discovered b y IRAS to be very brig ht at far-IR wavelengths, 
might also be not young and fairly nearby (IZuckerman et al.l 120071 ). In contrast, to the best 
of our knowledge, Spitzer has not yet discovered any nearby star not in a region of recent 
star formation nearly as infrared luminous as these four. By infrared luminosity we mean the 
fraction of a star's bolometric luminosity as seen from Earth that is absorbed and re-radiated 
by dust particles. For the above-mentioned four stars this fraction is in the range between 2 
and 75%. 

When the remarkable properties of BD+20 307 were first appreciated, a statistical anal- 
ysis of the frequency of occurrence of such very dusty stars, based on only one example , 
might be regarded as premature (and no such analysis was attempted by lSong et al.ll2005l ). 
However, with the recognition (below) that properties of HD 23514 are quite similar to 
those of BD+20 307, the phenomenon has been transformed from a miracle into a statistic. 
Consequently, following description of our observations, in Section 3 we discuss the occur- 
rence frequency of the very dusty phenomenon and what it might imply for the evolution of 
planetary systems in orbit around adolescent-age main sequence stars. 



2. Observations and Results 



The large mid-infrared excess of HD 23514 (Fig. 1) was discovered by IRAS only in 
its 12 iim band. Due to its large beam size, IRAS measurements often included many con- 
taminating background sou rces and some IRAS identified IR ex cess stars were subsequently 
found to be false positives (ISong et al.ll2002l : iRhee et al.ll2007bl ). In addition, all previously 
known dusty main sequence stars with excess emission detected at IRAS 12 yum, also had 
excess emission in at least one of the longer IRAS 25, 60, or 100 /im bands. This is because 
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the IR detection of cold Kuiper Belt analogs has been much more frequent than detection 
of warm asteroid belt analogs. Thus, significant ambiguity among IRAS excess candidates 
and lack of known main sequence stars with strong mid-IR excess emission perhaps helped 
to prolong the overlooking of the IRAS 12 /im detection of HD 23514. 



Spangler et al.l (120011 ) reported a marginal detection of dust excess emission from HD 23514 
at 60 and 90 /xm with a pointed observation of the Infrared Space Observatory (ISO). They 
did mention the IRAS 12 /im measurement, but then used only ISO data leading them to 
an incorrect conclusion about the dust properties of this star (Tdust ~ 70 K, L^ust/L^ ~3 
xlO~^). HD 23514 was rediscovered as a potential hot dust star from our search of main 
sequence mid-IR excess stars using pubhc Spitzer data (Rhee, Song & Zuckerman, in prepa- 
ration). Spitzer MIPS 24 /xm images have a fairly large field of view (5'x 5') and many 
field stars appear in most MIPS images. HD 23514 wa s included, serend ipitously, in the 



FEPS [Formation and Evolution of Planetary Systems, (IMeyer et al.l 120041 )] field of Pleiad 
HH 1182. 

Follow-up imaging observation of HD 23514 was carried out with the Near InfraRed 
Imager (NIRI) and Mid-IR Imager/Spectrometer (Michelle) at Gemini North Telescope. 
L' (3.8 /im) and M' (4.7/im) images were obtained with NIRI using a four-point dithered 
pattern. The standard "beam switching" mode was used for six mid-IR narrow-band images 
with Michelle by chopping the secondary at 2.7 Hz and nodding the telescope every ~30 
sec. For L' and M' images, dark frames were first subtracted from raw frames. After sky- 
subtraction, images were then flat-fielded using a sky frame made by median combination 
of dithered images. Images at each band were shifted, added, and averaged to produce a 
final image at each wavelength. For six mid-IR images, raw images were sky-subtracted 
using the sky frame from each chop pair. Subtraction of a nodded pair removed the thermal 
emission from the telescope. Standard stars, HD 22686 and HD 18884, were observed close 
in time and position to our target and used for absolute flux calibration of L' and M' images 
and of six mid-IR images, respectively. Finally aperture photometry was performed on both 
HD 23514 and the standard stars to compute flux density at each band. We used aperture 
radii of 0.96" and sky annuli of 1.42" & 1.97" for both target and the standard star. 

These ground-based images at 3.8 — 11.7/im bands show only one object at the expected 
target location (Gemini blind pointing accuracy is good to < 1"), thus verifying that the 
dust emission shown in Figs. 1 and 2 originates from HD 23514. Flux measurements at 
those wavelengths conflrmed excess emission above the stellar photosphere. 

Near-IR and mid-IR spectra of HD 23514 were o btained using the Ne ar-infrared Cross- 



dispersed Echelle Grating Spectrometer (NIRSPEC, iMcLean et al.l Il998l ) at Keck II Tele- 



scope and Michelle at Gemini North Telescope. NIRSPEC was used in a low resolution 
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mode (R ~ 2000) with the 42" x 0.570" sht to obtain a KL (2.8 - 3.7 //m) band spectrum 
of HD 23514. For N-band spectra, the low-resolution spectroscopic mode (R ~ 200) of 
Michelle was used with a 2 pixel-wide (0.402") slit. The N-band filter with a central wave- 
length of 10.5 /im was selected to give wavelength coverage of 7.7 — 14/im. REDSPEC, an 
IDL based reduction package for NIRSPEC, was used for the reduction of the NIRSPEC 
KL spectrum. Dark frames were first removed from raw frames. Using the sky frame of the 
nod pair for KL spectrum and of the chop pair for N-band spectrum, the resultant frames 
were sky-subtracted and flat-fielded. After the spectra of the standard stars (HD 210501 
for KL-band and HD 18884 for N-band spectra) were divided by Planck curves with each 
star's effective temperature (6400 K for HD 210501 and 3400 K for HD 18884), these ratioed 
spectra were then divided into the spectra of HD 23514 to remove telluric and instrumental 
signatures. Wavelength calibration was performed using Argon lamp spectrum for the NIR- 
SPEC KL spectrum and using atmospheric transition lines from an unchoppcd raw frame 
for Michelle N-band spectrum. Finally, photometry values at NIRI L' band and Michelle 
8.8 /xm band were used to fiux-calibrate KL and N-band spectra, respectively. Near-IR & 
mid-IR photometry of HD 23514 from both ground-based and space observations is listed in 
Table 1. 

We note that MSX detected HD 23514 in its A-band (8.28 /xm isophotal wavelength) 
with a catalog flux of 158 mJy (color uncorrected) significantly below our measured fluxes 
from the ground with Michelle narrow bands. Furthermore, the reported position of the 
MSX source is ~12" west of HD 23514. To reconcile the discrepancy in flux and check the 
nature of the ~3 sigma positional offset, we downloaded a 2x2 degee MSX A-band image of 
HD 23514. The image contains about two dozen MSX sources. When we overplotted MSX 
catalog positions on the MSX A-band image, only HD 23514 shows a mysterious shift of 
~12" from the obvious source in the image. All other MSX catalog positions fall right on 
bright sources in the image. Given the fact that no source other than HD 23514 appeared in 
the 32" X 24" field of Michelle images, we attribute the mysterious 12" offset to an erroneous 
astrometric correction of MSX. Although the nominal wavelength of the MSX A-band is 
8.28 /im, its effective wavelength is dependent on the source spectrum (hence the need for 
color correction). When a true source spectrum is very different from the assumed one 
{Fx ~ A~^) as in our case, the effective wavelength of the MSX A-band can be significantly 
shifted. To quote an MSX measured flux at its nominal wavelength of 8.28 fim, a color 
correction needs to be applied which can account for the apparent discrepancy between our 
narrow band and MSX catalog fluxes. 

In Fig. 1, the slope of the NIRSPEC KL spectrum agrees with the dust continuum flt we 
derive below. The 8 — 13 /xm spectrum. Figures 1 and 2, reveals warm small dust grains near 
HD 23514 through a prominent emission feature. For young stars and debris disks, the most 
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prominent spectral feature in the N-band is silicate emission. However, our N-band HD 23514 
spectrum peaking at ~9 fim is different from almost all other frequently seen silicate features 
that peak at 9—11 /im due to various combinations of olivine, pyroxene, and other minerals 
(Fig. 2). Possible carriers of this bizarre 9 fim emission feature among common minerals 
in our Solar System and Earth's surface are tektosilicates and sulfates. Tektosilicates are 
a group of light colored silicate minerals and this group contains most common minerals 
(quartz, feldspar, etc.) seen on Earth's surface. About 75% of Earth's crust is composed of 
tektosilicates. However, explaining the strong 9 /im feature without accompanying prominent 
olivine and pyroxene signatures is challenging. For example, it is difficult to imagine an 
extraordinary amount of sulfates at HD 23514 over more commonly appearing minerals 
unless the chemical composition of HD 23514 is very different from solar. The huge quantity 
of dust needed to match the HD 23514 SED can be generated by catastrophic collisions 
among planetary embryos or even a planet-planet collision (see § 3). The latter mimics 
the postulated Moon-creating collision between th e young Earth and a Mars size planet 
( iHartmann fc David Il975l : ICameron fc Ward! Il976l ). Crustal material ejected from such a 
collision may naturally explain our HD 23514 N-band spectrum. However, it is hard to 
explain how crustal material was ejected predominantly over mantle material considering 
that Earth's crust is a thin layer occupying only <5% of volume compared to the mantle. 
Nonetheless, the unusual N-band spectrum of HD 23514 must bear a clue to the origin of dust 
and a wider range mid-IR spectrum is needed for more detailed analysis. At this stage, we 
likely rule out the case of collisional grinding of many asteroids as the source of dust around 
HD 23514 since olivine and pyroxene should be the dominant mine rals in such environm ents. 
We note that few objects out of 111 T Tauri stars investigated by lFurlan et al.l (120061 ) . who 
used SpitzerlRS, show a mid-IR emission feature peaking near 9 /im (e.g. IRAS 04187-1-1927 
and CZ Tau), as found in HD 23514. 



3. Discussion 

Excess emission peaking at mid-IR wavelengths indicates that dust must be warm and 
close to the central star. The temperature, as well as the amount of dust and its distance from 
the central star, is constrained by creating a Spectral Energy Distribution (SED) assuming 
that dust exists as an optically thin ring. We produced an SED of HD 23514 (Fig. 1) by 
f itting observed rneasur ements at optical and infrared bands with a stellar photosphere model 
faauschildt et al.lll999h and a single temperature blackbody of T = 750 K. Large blackbody 



grains in thermal equilibrium at 750 K would be located ~0.25AU from HD 23514. Even 
small grains that radiate less efficiently, especially those responsible for the mid-IR em ission 



feature, likely lie within a few AU of the central star. We noted in iRhee et al.l (l2007al ) that 
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stars with warm dust emission, ages between 10 and 30 Myr, and spectral types from GO 
to A do not show any evidence in their SEDs of the presence also of cold dust. In such 
cases, even the small dust particles that carry the strong mid-IR emission feature are likely 
located close to the stars. Likewise, there is no obvious evidence for cold dust associated 



with BD+20 307 (iSong et al.l l2005l : A. Weinberger et al. 2007, in preparation). However, 
this may not be the case for HD 23514 because Spitzer and ISO points in its SED all lie 
somewhat above the 750 K dust continuum as shown in Fig. 1, thus suggesting the presence 
of cooler dust further from the star. Still, caution is appropriate in interpretation of the 
Spitzer 25 /xm flux density measurement as it may be elevated by inclusion of the red wing 
of a silicate emission feature (e.g., A. Weinberger et al. 2007, in preparation). And the 60 
and 90 fim ISO points, respectively, lie only 2 and 2.3 sigma above the 750 K dust continuum 
line in Fig. 1. A Spitzer 70 fim flux measurement is highly desirable to clarify the presence 
of cool dust. 

A standard method for characterizing the amount of dust orbiting a star is through the 
quantity r (= L//j/L*) where L//j is the excess luminosity above the photosphere emitted 
at infrared wavelengths and L,, is the bolometric luminosity of the star. We obtained r ~ 
2 X 10~^ by dividing the infrared excess between 2.3 /im and 90 /zm by the stellar bolometric 
luminosity (2.8 Lq); this is ~ 10^ times greater than that of the Sun's current zodiacal cloud 
(r ~ 10~^). HD 23514 thus joins BD+20 307 as the two Sun-like main sequence stars with by 
far the largest known fractional ir ifrared luminosit ies (Table 2). The age of BD+20 307 is at 



least a few hundred million years (ISong et al.ll2005l ). As a member of the Pleiades (HII 1132) 



the age of HD 23514 is ~100Myr. There can be little doubt of cluster membership because 
in the plane of the sky HD 23514 is located well inside the clustei0, sharing common proper 
motion, and its radial velocity of 5. 9+0.5 km s~^ is in good agreement with the velocity of 



the Pleiades cluster, 6.0+1.0 km s ^ (ILiu et al.l 1 19911 ). 



Currently only a handful of stars with ages > 50 Myr show warm excess emission 
(T ^ 150K), indicative of planetesimals in the terrestrial planet zone. Table 2 lists some 
parameters of these stars. BD+20 307 and HD 23514 stand out among them with very high 
dust temperature (T ^ 600 K) and r > 10~^. The remaining four warm excess stars have 
cooler dust and r ~ 10~^. 

The SEDs of HD 23514 and BD+20 307 exhibit excess near-IR emission beginning at 
wavelengths ~4 fim. Given the potential importance of stellar wind drag on the dust particles 
(see below), the absence of really hot dust (>1000K), suggests the possible presence of a 



picture of the Pleiades Cluster with HD 23514 indicated is available from the Gemini Observatory 
Web site, http://www.gemini.edu/index.php?option=content^ task=viewfcid=259. 
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"sweeper planet" closer to the stars than the dust. Such a situation pertains at HD 69830 
where, also, no very hot dust is seen and where Neptune-mass plane ts interior to the dust 
disk are known to exist from precision radial velocity measurements ( iLovis et al.ll2006l ). 



Using a flat disk model. lJural (120031 ) and I Jura et al.l (120071 ) have successfully reproduced 
the IR emission from flat, geometrically thin, dust disks orbiting some white dwarfs. This 
geometry implies the absence of signiflcant gravitational perturbations by objects with sub- 
stantial mass located in the vicinity of the dust. In contrast, we flnd that a flat disk of dust 
particles generates insufficient mid-IR flux to match the SEDs of HD 23514 and BD+20 307. 
Thus, the dust orbiting these stars is puffed up in the vertical direction, perhaps as a result 
of the gravitational field of the above mentioned sweeper planet, or the gravity of planetary 
embryos as discussed below, or both. 

Given the ages of HD 23514 and BD-l-20 307, it is natural to ask whether their huge 
warm dust burdens were generated by events analogous to those that occurr ed during the 



"Late Heavy Bombardment" (LHB) in our solar system. One current model (IGomes et al. 



20051 ) attributes the LHB to a rapid migration of the giant planets that destabilized the 
orbits of objects in the Kuiper Belt and the rnain a steroid belt hundreds of millions of years 
after the formation of the Sun. IWyatt et al.l (120071 ) have proposed a similar model for most 
of the stars listed in our Table 2, with the difference that they strongly favor an origin of 
the parent bodies in a region more analogous to the cold Kuiper Belt than the asteroid belt. 
We note, however, there is no evidence for cold dust at BD-(-20 307, ( Lep, and HD 69830, 
cold dust that might reasonably be expected at stars with so much warm dust should all 
parent bodies originate in distant cold regions. The situation is more ambiguous for a Table 
2 star like r] Crv with c lear e vidence for substantial amounts of cold dust (see Fig. 3). 
Furthermore, ICuk et al.l (120061 ) argue that LHB was a localized Earth- Moon system activity 
rather than a global, inner solar system, event. Thus, because the cause (s) of the LHB 
remain unsettled, we do not further pursue a relationship between LHB and the high r 
warm excess phenomenon. 

Initially, when there was only one known main sequence star with very large r (> 10^^), 
the BD+20 307 phenomenon might have been regarded as a "miracle" , so that a statistical 
analysis of the occurrence rate would have been of questionable value. Now, however, with 
HD 23514, the frequency of occurrence of such extraordinarily dusty stars can be treated 
statistically more reliably. BD+2 307 is a field st ar with estimated age similar to that of 
the Ursa Ma j oris movin g group (jSong et al.l 120051 ). whose age is probably about 400 Myr 
(IZuckerman et al.ll2006[ ). IRAS was sufficiently sensitive to detect main sequence G-type 



stars with r > 10 ^ out to ~150 pc. There are ~18400 Hipparcos dwarfs with spectral 
types between F4 and KO out to 130 pc, the distance to the Pleiades. But most of these 
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are old stars. If the star formation rate was approximately uniform during the past 5 Gyr 
years, then there are about ~1800 Hipparcos dwarfs of age ~500 Myr out to 130 pc. Then 
a Hipparcos dwarf with r > 10"^ is found among solar-type field stars about one time out 
of 1500 (after dropping ~300 Hipparcos-measured members of nearby stellar clusters). 

In the solar vicinity several stellar clusters have solar-type members. We select four rich 
nearby clusters with ages 70—700 Myr: Hyades, Pleiades, a Persei, and Praesepe. In these 
stellar clusters, there are about ~400 dwarfs stars with spectral types between F4 & KO 



dde Bruiine et al.ll200ll : IStauffer et al.ll2007l : iLodieu et al.ll2005l : lAdams et al.ll2002l ). Among 
these only HD 23514 is identified with r > 10~^. Therefore the occurrence rate of a dwarf 
with r > 10~^ in the nearby stellar clusters is, at most, about one out of 400. Combining 
this result with that for BD+20 307, indicates that the very high r, warm dust, phenomenon 
manifests itself in about one adolescent star (age a few 100 Myr) in 1000. If all F4-K0 stars 
display this phenomenon as adolescents, then the lifetime of the phenomenon at a typical 
solar-like star is a few 100,000 years. 

To interpret our observations, we consider a mo del of colliding planetary embryos. In 
a ser i es of papers, Agnor, Asphaug, and colleagues (lAgnor et al.l Il999l : lAgnor fc Asphaug 
2004J : lAsphaug et al.ll2006l ) considered the collisions of large bodies in the late stages of the 
formation of planets in the terrestrial planet zone. Based on their models and those they 
attribute to earlier researchers (e.g., G.W. Wetherill), we may draw the following conclusions. 
The process of terrestrial planet formation involves the formation of a minimum of many 
hundreds of planetary embryos of dimensions >1000 kilometers. These collide and either 
they coalesce or, oftentimes, the smaller embryo fragments into smaller objects along with 
the ejection of "copious debris". While the mass spectrum of the fragments is not well 
constrained, no large monoliths survive following disruption of large solid bodies. Rather a 
typical large fragment size might be ~100m. Collisions of planetary embryos continue for 
as long as a few 100 Myr, i.e., to the ages of HD 23514 and BD+20 307; in the following 
discussion, we assume these to be 100 and 400 Myr, respectively. Additional discussion of 
catastrophic fra gmer itation of planetary s ystern . bodies of moderate size may be found in 
Fuiiwaral Jl98oh and iHousen fc Holsapplel fll99oh . 



For these assumed ages and an occurrence rate of one in 1000 stars, we will use a 
lifetime of 250,000 years for the HD 23514/BD+20 307 phenomenon at a typical adolescent- 
age solar-type star. Small particles now in orbit around these two stars will be lost in a much 
shorter time span and must be replenished many times over. One possible loss mechanism 
is a collisional cascade that breaks particles down in size until, when their radii become 
as small as a few tenths of a micron, they become subject to radiation pressure blowout. 
Other loss mechanisms are Poynting-Robertson (PR) and stellar wind drag. As mentioned 
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above, we assume that the initial mass spectrum is a result of the collision of two planetary 
embryos, but the spectrum of the collision fragments is not well characterized. Therefore, 
we assume that collisions are sufficiently frequent to establish an approximately equilibrium 
size distribution: 

N{a)da = Noa'^'^da (1) 
where N(a) is the number of particles per cm^ with radii bet ween a and a + da (see, e.g.. 



Dohnanyilll969l : IWilliams fc Wetherilllll994l : IChen &: Jurall200ll ). The smallest particle radius 



in this distribution may be set by radiation pressure blow-out; for the mass (1.35 M©) and 
luminosity (2.8 L©) we estimate for HD 23514, this radius is ~0.5 microns. 

With this size distribution most of the mass (M) is carried by the largest particles while 
most of the surface area (r) is due to small particles with radii not much larger than the 
submicron-size blowout radius. Specifically, 

Mr. [ ^a'a--'-'da ~ a'^' (2) 



3 

r(a) ~ / N{a)'Ka^da ~ / Tra^a'^'^da ~ a~^^^ (3) 



and 

= P/t ~ a'/' (4) 

where tc is the coUisional lifetime and P is the orbital period. Thus, the smallest particles 
collide the fastest, for HD 23514 in about 50 years at lAU. Larger objects take longer to 
collide destructively. They will then be broken down into smaller fragments such that after a 
collision of two roughly equal mass objects the largest left-over fragment has a radius about 
1/2 that of a collider (S. Kenyon 2007, private comm.). As more mass is carried by the 
larger colliders, with the above N(a), the ratio of mass to collision time is independent of a, 
and there will thus be a supply of material approximately constant with time as the largest 
objects are eventually whittled down to micron size dust. In the Agnor/Asphaug picture 
outlined above, the largest initial fragments of a collision of planetary embryos might have 
a ~ 100 m. Thus, if the lifetime of a = 1 fim particles is ~50 years due either to collisions 
or to stellar wind drag (see below), then the lifetime of 100 m fragments will be ^500,000 
years. Thus, catastrophic disruption of large planetary embryos (see mass estimate below) 
can supply material for a time equal to the 250,000-year event lifetime indicated by our 
observations. 

In addition to collisions, PR and stellar wind drag are loss mechanisms for small dust 
part icles. The tirn e scale for PR drag at HD 23514 may be evaluated using, for example, Equ. 



5 in lChen fc Jural (120011 ). We assume L,,, = 2.8 L© (corresponding to a distance of 130 pc to 
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HD 23514), a dust particle orbital semi- major axis of one AU, and a radius and density of a 
typical individual grain to be 1 fim and 2.5 g/cm^, respectively. Then the lifetime against PR 
drag is ~1000 years. For the present day Sun, PR drag acts ab out 3 times naore r apidly than 
drag due to the solar wind (jPlavchan et al.ll2005l ). However, IWood et al.l (120021 ) estimated 
that the winds of solar type stars decline as time to the 2.00 ± 0.52 power beginning at 
ages ~10% of the current age of the Sun. More recent observation and analys is somewhat 



cloud quantitatiy e representation of the wind strength as a function of time (IWood et al. 



20051 : IWoodll2006l ). Based on these references, we assume that, for adolescent-age stars like 
HD 23514 and BD+20 307, wind drag will dominate PR drag by a factor 10 - 30. Then the 
lifetime of the small particles considered above against wind drag will be only ~50 years, 
i.e., comparable to the collision times. Because PR and wind drag times are proportional 
to a, while the collision time is proportional to a^/^, the orbital lifetimes of large particles 
(rocks) are determined by collisions. 

To determine how rapidly mass is lost due to either stellar wind drag or collisions we 
estimate the minimum dust mass (M^m) needed to intercept 2% of the light emitted by 
HD 23514. This is 

Mmin = 167rrapi?V3 (5) 

Where r = 0.02, a = 1 fxm, p = 2.5g/cm^ and R = lAU. Then M^in ~ 2 x 10^^ g, with 
a corresponding mass loss rate, M ~ 10^'^ g/s. For BD+20 307, r is twice as large, so M 
is ~2 xlO^^ g/s. In 250,000 years, the total mass lost per star will be ~10^^ g— the mass 
of Earth's moon. Or, for the above assumed average density, an object with radius ~2000 
km. Of course, this mass need not all be produced in one single catastrophic collision, but 
might rather be a consequence of multiple collisions of smaller planetary embryos spaced 
over hundreds of millions of years. 



4. Conclusions 



We show that substantial quantities of warm, small, dust particles orbit HD 23514, a 
solar-type mem ber of the 100 M yr old Pleiades cluster. A similar phenomenon was previ- 
ously reported (ISong et al.ll2005l ) for the somewhat older, solar-type, field star BD+20 307. 
Models for catastrophic collisions of planetar y embryos orbiting ~100 Myr old stars (e.g., 
Kenyon fc Bromleyll2005l : lAsphaug et al.ll2006l ) will naturally produce such warm dusty disks. 
Our data are consistent with these model predictions, provided that such catastrophic events 
followed by a subsequent coUisional cascade convert of order the mass equivalent of Earth's 
Moon to tiny dust particles during the early lifetime of many (perhaps most) Sun-like stars. 
For example, in the case of our solar system, by itself the collision that is postulated to have 
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generated the Moon likely would have sent a comparable mass of debris into interplanetary 
orbits. Infrared data for stars such as HD 23514 and BD+20 307 are consistent with and 
may well validate the standard picture of violent formation of terrestrial-like planets in the 
early years of planetary systems. 
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Table 1. HD 23514 near-IR & mid-IR Photometry 



Filter 


Central Wavelegth 


Flux Density 


Uncertainty 






(/im) 


(mJy) 


(mJy) 


Instrument 


B 


0.44 


466^ 


llO*" 


TYCHO-2 


V 


0.55 


600^^ 


86^ 


TYCHO-2 


T 

J 


1.25 


670 


15 


2MASS 


TT 

H 


1.65 


516 


11 


2MASS 


K, 


2.20 


369 


9 


2MASS 


L' 


3.78 


196 


10 


Gemini NIRI 


M' 


4.68 


175 


7 


Gemini NIRI 


Si-1 


7.7 


197 


10 


Gemini Michelle 


Si-2 


8.8 


234 


12 


Gemini Michelle 


Si-3 


9.7 


209 


10 


Gemini Michelle 


Si-4 


10.3 


140 


7 


Gemini Michelle 


Si-5 


11.6 


104 


5 


Gemini Michelle 


Si-6 


12.5 


95 


5 


Gemini Michelle 


12/im 


11.5 


184 


26 


IRAS 


24/im 


24.0 


66.5 


2.7 


MIPS 


60/im 


60.0 


28 


10 


ISO 


90/um 


90.0 


26 


10 


ISO 



^The standard B and V magnitu des were obtai ned by converting Tycho B 
and V magnitudes using Table 2 in iBesselll (120001 ). 

^B- & V- flux density uncertainties were computed assuming 0.200 for their 
magnitude uncertainties in order to compensate for some missing opacity 
species in the model spectrum (see lRhee et al.ll2007bl .) 
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Table 2. Main sequence stars with debris systems in the terrestrial planetary zone 



Dust 

Spectral Temperature r Age Cold 

Object Type (K) (xlO~^) (Myr) dust References'' 



HD 23514'' 


F6V 


750 


200 


100 


Maybe 


1 


BD+20 307^ 


GOV 


650 


400 


400 


no 


2 


C Lep 


A3 


190 


0.65 


300 


no 


3,4 


HD 72905'^ 


G1.5 


? 


1 


400 


yes 


5,6 


T] Corvi 


F2V 


180 & 30 


5 


600 


yes 


4,6,7 


HD 69830^^ 


KOV 


? 


2 


2000 


no 


6,8 



Note. — We define the terrestrial planet zone (TPZ) to be the region where 
dust parti cles that radiate l ike blackbodies will attain a temperature of at least 
150 K fsee lRhee et "al]l2007al b 



mm 5. 



Beichm 



Beichman et al 



an et al. 



1. T his paper 2. ISong et al.l (120051) 3. IChen k Jural (l200ll'l 4 



((20061) 6. IWvatt et al.l ([20071) 7. IWvatt et al 



Chen et al. 
(|200,5l ) 8 



(120051) 

''Our fit to the SED in Fig. 1 implies a photospheric temperature of 6400 K 
and a stellar radius = 1.28 i?© for an assumed distance to HD 23514 of 130 pc. 
However, based on pre-main sequence evolution models of Baraffe et al. (1998,2002), 
a ~100Myr old 6400 K star has radius 1.38 i?© and mass 1.35 M©. If HD 23514 
is a single star, then the discrepancy between the two radii would be eliminated 
if the actual di stance to HP 235 14 is ^ 140 pc. Th e spectral type of HD 23514 is 
listed as F5 in ICrav et al.l ([20011) while \Cod (|2000l ) gives F7 for a main sequence 
star with T = 6400 K. We adopt F6 in this paper. 



=T = 6000 K, R = 1.25 i?0 (jSong et al.ll2005l ) 



'^Following IWvatt et al.1 (120071 ) we list HD 69830 k HD 72905 as potenti ally 
having dust particles in the TPZ. lBeichman et al.l (|2005l ) and lLisse et al.l (|2007l ) fit 
a complex model to the mid-IR spectrum of HD 69830 and derive an underlying 
dust continuum temperature of 400 K. However, given the number of free param- 
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ctcrs included, and not included (e.g. particle shape), in these models we regard 
the dust temperature as not well constrained. 



-17- 



B V J H K L' M' 12|j,m 25|J,m 60|J,m90|J,m 




1 10 100 

Wavelength (jim) 



Fig. 1. — Spectral Energy Distribution of HD 23514. NIRSPEC spectrum is shown in yellow 
(wavelengths 2.8 — 3.7/im), while the Michelle spectrum is in maroon (7.8 — 13.3/xm). The 
derived stellar parameters for HD 23514 are given in a note to Table 1. For each measurement 
the horizontal bars indicate the passband of the filter used and the vertical bars depict 1 a 
flux uncertainties. 
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Fig. 2. — Michelle mid-IR spectrum of HD 23514. A 2 " wide sky region ~3 " from the HD 
23514 spectrum was used to compute the uncertainty at each wavelength. Representatve 1 a 
uncertainties of the spectrum (vertical bars) are given at each half micron. For c omparison. 



the m id-IR, dotted-line, spectrum of BD+20 307 is reproduced from Fig. 1 in ISong et al. 



( I2OO5I ). The spectrum of HD 23514 has almost twice the spectral resolution as the Keck LWS 
spectrum of BD+20 307. For each photometry measurement (green square) the horizontal 
bars indicate the passband of the filter used and the vertical bars depict 1 a uncertainties. 
As noted in § 2, the shape and location of the BD+20 307 feature is similar to that unsually 
seen in pre-main sequence stars and solar system comets. 
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B V 



J H K 



12 



25 




1 — I — I — I — r~r 



60 



450 850 



"T — I — I — r~r 



T 1 — TTTT 



• JHK{2MASS) 

■ 'RAs ; 

A SpitzerMIPS 24 & 70 |im 

— Spitzer ;IRS spectrum 

• SGUBA 450 & 850 |im 

— 6800K Stellar photosph; 




star + dust continuum (180 K + 30K) 



J I I I -I I 
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J I I I'll 



10 



100 



1000 



Wavelength (jim) 



Fig. 3. — Sp e ctral energy distribution of r] Crv. This should replace the SED in Fig. 4 in 
Wyatt et al.l (120051 ) because: (1) we plot IRAS FSC data in preference to the less accurate 



IRAS PSC data used by Wyatt et al (2005), and (2) we plot the Spitzer 24 micron MIPS 
measurement. Note that, as a result of these changes, our decomposition of the infrared 
portion of the SED is very different from that of Wyatt et al. Christine Chen kindly provided 
the fiux-calibrated Spitzer IRS spectrum of rj Crv. 



